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The new version (INCL4.6) of the Liege intranuclear cascade (INC) model for the description of 
spallation reactions is presented in detail. Compared to the standard version (INCL4.2), it incorpo- 
rates several new features, the most important of which are: (i) the inclusion of cluster production 
through a dynamical phase space coalescence model, (ii) the Coulomb deflection for entering and 
outgoing charged particles, (Hi) the improvement of the treatment of Pauli blocking and of soft col- 
hsions, (iv) the introduction of experimental threshold values for the emission of particles, (v) the 
improvement of pion dynamics, (vi) & detailed procedure for the treatment of light-cluster induced 
reactions taking care of the effects of binding energy of the nucleons inside the incident cluster and 
of the possible fusion reaction at low energy. Performances of the new model concerning nucleon- 
induced reactions are illustrated by a comparison with experimental data covering total reaction 
cross-sections, neutron, proton, pion and composite double differential cross-sections, neutron multi- 
plicities, residue mass and charge distributions, and residue recoil velocity distributions. Whenever 
necessary, the INCL4.6 model is coupled to the ABLA07 de-excitation model and the respective 
merits of the two models are then tentatively disentangled. Good agreement is generally obtained 
in the 200 MeV-2 GeV range. Below 200 MeV and down to a few tens of MeV, the total reaction 
cross section is well reproduced and differential cross sections are reasonably well described. The 
model is also tested for light-ion induced reactions at low energy, below 100 MeV incident energy 
per nucleon. Beyond presenting the update of the INCL4.2 model, attention has been paid to ap- 
plications of the new model to three topics for which some particular aspects are discussed for the 
first time. The flrst topic is the production of clusters heavier than alpha particle. It is shown that 
the energy spectra of these produced clusters are consistent with coalescence. The second topic 
regards the longitudinal residue recoil velocity and its fluctuations. Excellent results are obtained 
for these quantities. It addition, it is shown that the distributions of these quantities display typical 
random- walk characterics, at least for not too large mass losses. They are interpreted as a direct 
consequence of the independence of successive binary collisions occurring during the cascade process. 
The last topic concerns the total reaction cross section and the residue production cross sections 
for low energy incident light ions. It is shown that our new model can give a rather satisfactory 
account of these cross sections, offering so an alternative to fusion models and the advantage of a 
single model for the progressive change from fusion to pre-equilibrium mechanisms. 



PACS numbers: 25.40.-h,25.40.Sc,25.45.-z,25.55.-e 



I. INTRODUCTION 



It is largely accepted that nucleon-induced spallation 
reactions proceed through a two-stage process: a first 
stage dominated by hard nucleon-nucleon (NN) collisions 
emitting fast particles, followed by the de-excitation of 
a more or less thermalized remnant, akin to evapora- 
tion and/or fission. The commonly used tools to de- 
scribe these reactions result from the coupling of an intra- 
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nuclear cascade (INC) model for the first stage to an 
evaporation-fission model for the second stage. There 
exist many INC and de-excitation models, differing from 
each other by the ingredients and by the results. In the 
last years, there has been a very strong development 
of applications involving spallation reactions, including 
neutron sources for condensed matter and material stud- 
ies [H, transmutation of nuclear waste simulation of 
experimental set-ups in nuclear and particle physics 
production of rare isotopes [5,], protection against radi- 
ation near accelerators and in space missions Q, inter- 
action of cosmic rays in the atmosphere Q and cancer 
hadrontherapy [1]. This activity is at the origin of a 
strong demand for the improvement of the accuracy of 
the reaction models and for reliable data to benchmark 
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these models efSciently. Indeed, in Europe, an important 
effort lias been devoted in the last fifteen years to col- 
lect and produce high-quality data concerning all emis- 
sion channels: neutrons, light charged particles (Icp) and 
residues d, [l^l- Meanwhile, model developers have con- 
stantly attempted to improve their models. We are in- 
terested here in the development of the Liege INC model 
(INCL4). The first standard version (INCL4.2) is de- 
scribed extensively in Ref. Illl . Coupled with the ABLA 
evaporation/fission model jl3 - (l4| (actually the KHSvSp 
version), it is able to describe fairly well a large body of 
experimental data in the 200 MeV-2 GeV range, namely 
total reaction cross sections, double differential proton, 
neutron and pion production cross sections, residue mass 
and charge spectra, isotopic distributions and, to some 
extent, recoil energies [111] . This was achieved without 
any parameter adjustment in the INC model and with the 
use of standard values for the parameters of the ABLA 
model. More precisely, there is no fitting parameter in 
the INCL4.2 model. The basic nuclear physics param- 
eters, like the radius and the diffuseness of nuclei, are 
taken from phenomenology and "technical parameters" , 
like those used to evaluate phase space occupation enter- 
ing the Pauli blockers, have been fixed once for all. In ad- 
dition the stopping time is generated self-consistently by 
the model itself. Nevertheless, the results obtained with 
INCL4.2 show some systematic shortcomings and/or dis- 
crepancies with experimental data. They have been iden- 
tified in the frame of the European HINDAS and NUDA- 
TRA collaborations [1, ^ and are shortly presented be- 
low. The most important of these shortcomings is the 
inability of INCL4.2 to produce Icp's (except protons of 
course). Since the release of the INCL4.2 model, a con- 
stant effort has been made to palliate its deficiencies. 

It is the purpose of this paper to set out the present sta- 
tus of the INCL4 model, that is designated as INCL4.6. 
The various points of improvements are described be- 
low. But, this paper is not a mere update of the INCL 
model. It highlights at least three innovative points con- 
cerning spallation reactions and bearing on (i) produc- 
tion of clusters in the cascade, (ii) residue recoil veloci- 
ties and (iii) the behaviour of the model at low energy, 
respectively. Concerning the first point, a dynamical 
coalescence model has been implemented to allow the 
emission of clusters in the cascade stage. An interme- 
diate version of our model [l5| . sometimes denoted as 
INCL4.3, consisting basically of INCL4.2 plus the treat- 
ment of the emission of light clusters, up to alpha parti- 
cles, has opened the path to the improvement along this 
line. We report here on the continuation of this effort 
and we show, for the first time, that emission of heav- 
ier clusters (in practice, up to A=8) can be explained by 
dynamical coalescence. Concerning point (ii), we show 
below that we are now able to reproduce average values 
and standard deviations of the residue velocity distribu- 
tions. Finally, let us come to point (iii). It is generally 
stated that INC models cannot be reliable below ^-^200 
MeV incident energy, although in some cases where such 



models are used, they occasionally give surprisingly sat- 
isfactory results. In the last years, we have tentatively 
but systematically improved the model in this energy 
range. We considerably improved the model concern- 
ing reactions induced by nucleons and by light clusters 
of nucleons, typically from deuterons to alpha particles, 
and this for incident energies per nucleon down to a few 
MeV. Our model is now able to produce good total re- 
action cross sections for both incident nucleons and in- 
cident light clusters. In addition, we show below that 
we are reproducing reasonably well residue production 
cross sections for cluster-induced reactions at low energy. 
The developments concerning the above-mentioned three 
points were motivated in part by studies of thick targets, 
concerning release of volatile elements, such as H and He 
isotopes (point (z)), radiation damages (point (ii)) and 
radiotoxicity (point {iii)). Sometimes, the latter may 
crucially depend upon secondary reactions induced by 
light clusters produced in primary collisions [16, 16]. 

The predictive power of this improved model is con- 
siderably better than the one of INCL4.2. This has been 
verified on an intermediate version (named INCL4.5, not 
very different from INCL4.6, and defined below) in the 
course of an intercomparison of spallation codes orga- 
nized by the IAEA [18I - I2Q ] . As one may surmise, the pas- 
sage from INCL4.2 to INCL4.6 (even to INCL4.5) could 
not have been done entirely in the same spirit of the 
building of INCL4.2, i. e. resting on known phenomel- 
ogy with well determined parameters. We were forced to 
introduce less well established features with less solidly 
determined parameters. The purpose of this paper is to 
present all the new features that have been introduced 
in INCL4.6 and to give, whenever possible, the physics 
motivation for these new features. Alth oug h a lot of re- 
sults with INCL4.5 are available in Ref. [I^, perhaps not 
sufficiently commented or detailed, we present here the 
most important ones, especially to illustrate the physics 
aspects of the improvements brought into the model. 

The paper is divided as follows. In Section II, we de- 
scribe the INCL model with emphasis on the new features 
of INCL4.6. Section III is devoted to an extensive com- 
parison with a representative panel of experimental data, 
for nucleon-induced reactions. In Section IV, we discuss 
theoretical results for light cluster-induced reactions con- 
cerning a few key experiments. Finally, in Section V, we 
critically examine the new features and their effects and 
we present our conclusion. 



II. DESCRIPTION OF THE INCL4.6 MODEL 
A. A short reminder of the INCL4.2 model 

The INCL4.2 model is extensively described in 
Ref. [ll|. It is sufficient here to remind the salient fea- 
tures. The INCL4 model is a time-like intranuclear cas- 
cade model. In the initial state, all nucleons are prepared 
in phase space. Target nucleons are given position and 
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momentum at random in agreement with a Saxon- Woods 
and a Fermi sphere distributions, respectively. They are 
moving in a constant potential well, the same for protons 
and neutrons, describing the nuclear mean field. The in- 
cident particle (nucleon or pion) is given the appropriate 
energy and an impact parameter at random. All nucleons 
are then set into motion and followed in space-time. The 
cascade process involves binary collisions between nucle- 
ons and produced pions and Delta resonances. Particles 
can be transmitted through, or reflected on, the surface 
of the square well potential they feel. Delta resonances 
decay according to their lifetime. Details on collisions, 
utilized cross sections, Pauli blocking, etc, may be found 
in Ref. [llj and previous publications cited therein. We 
remind that the stopping time of the cascade is deter- 
mined self-consistently by the model itself. It can simply 
be parametrized (in fm/c) by 

tstop = 29.8^^■l^ (1) 

for incident nucleons(we denote by Zt and At the charge 
and mass numbers of the target). 

It is useful to remind some other specific features which 
may be of importance for the rest of the paper. At the 
beginning of the cascade process, the incident nucleon 
or pion is located with its own impact parameter on the 
surface of the "working sphere" , centered on the target 
with a radius 

Rmax — Ra + 8a, (2) 

where i?o and a are respectively the radius and the dif- 
fuseness of the target nucleus density. Particles are mov- 
ing along straight-line trajectories between collisions in- 
side the working sphere. They are divided into partic- 
ipants and spectators in the usual sense. When partic- 
ipants leave the working sphere, they are considered as 
ejectiles and do not interact any more. Nucleons are mov- 
ing in a potential well with constant depth and with a 
radius which is dependent upon their momentum. The 
potential radius for particles with energy larger than the 
Fermi energy is also taken to be equal to Rmax- Pions 
do not experience any potential. Motivation and details 
are given in Ref. [ll| . 

The INCL4.2 model can accomodate light clusters (up 
to a-particles) as projectiles. In this case, the nucle- 
ons are given initially position and momentum inside the 
cluster and the latter is positioned at the beginning in 
such a way that one of its nucleons is touching the work- 
ing sphere. See Ref. [ll| for detail. 

We comment a little bit on soft NN collisions. The lat- 
ter do not contribute very much to the ejection of parti- 
cles in the cascade. There are several arguments based on 
nuclear transport theory indicating that these collisions 
should be disregarded, since they induce slow modifica- 
tions of the particle density distribution, a feature that 
is supposedly taken into account by the average nuclear 
potential (2l[-[25|. In INCL4.2, soft NN collisions, with 
cm. energy smaller than ^^=1925 MeV, are simply dis- 
regarded. This cut-off value may seem to be rather large. 



but one has to realize that low energy NN collisions are 
largely Pauli-blocked in spallation reactions. Lowering 
the value of ^/sq does not significantly change the results 
at high incident energy. It is no longer the case at low 
energy, say below ~200 MeV, as discussed in Ref. [26j . 
This matter will be re-examined in this paper. 



B. Shortcomings of the INCL4.2 model 

The shortcomings of the INCL4.2 model have been 
identified in various places [9l-[lT| . We just remind them 
very briefly. Some phenomenological aspects of nuclear 
physics are neglected. The model cannot accommodate 
production of clusters in the cascade, i.e. with a kinetic 
energy definitely larger than the typical evaporation en- 
ergies, as it can be seen experimentally. Concerning the 
predictive power of the model, several deficiencies can 
be noted. Pion production is generally overestimated. 
Quasi-elastic peaks in (p, n) reactions are generally too 
narrow and sometimes underestimated. Finally, reaction 
cross sections are severely underpredicted below ^100 
MeV. Residue production cross sections are sometimes 
unsatisfactorily reproduced, especially for residues close 
to the target. For matter of convenience, we separate 
below the new features of the INCL4.6 model into those 
which are included in the intermediate INCL4.5 version 
and those which are posterior to this version. 



C. Main new features in the INCL4.5 model 

1. Introduction of known phenomenology 

A. Isospin and energy- dependent potential well for the 
nucleons. The depth of the potential well felt by the 
nucleons is dependent upon the energy of the nucleons 
and is not the same for protons and neutrons. The en- 
ergy dependence is taken from the phenomenology of the 
real part of the optical-model potential [27l - [29| . Roughly 
speaking, the potential depth decreases regularly with in- 
creasing energy, from ordinary values at the Fermi level 
to zero at roughly 200 MeV. The isospin dependence is 
such that the neutron and proton Fermi levels have the 
same energy. For more detail, see Ref. (s^- The influ- 
ence of this modification is relatively small, except for 
special quantities, like the production cross sections for 
isotopes with an extra unit of charge compared to the 
target [30|.[3H. 

B. Average potential for pions. An average isospin- 
dependent potential well, of the Lane type [32| , is intro- 
duced for pions, as well as reflection and/or transmission 
at the border of this potential. The depth of the potential 
has been taken, as far as possible, from the phenomenol- 
ogy of the real part of the pion-nucleus optical potential 
(dispersive effects due to the strong imaginary part have 
to be removed). This depth amounts to 22 MeV for 7r''"'s 
and 38 MeV for 7r~'s on a Pb target. The radius of the 
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potential is taken as Rq + 4a, in rough accordance with 
phenomenology. This modification and its effect are pre- 
sented in detail in Ref. [13] ■ In general, it reduces the 
pion production cross section, mitigating so the overesti- 
mate by INCL4.2, as iUustrated below by Fig. [TTl 

C. Deflection of charged particles in the Coulomb field. 
Once an impact parameter is selected for the incident nu- 
cleon, the cascade process is initiated with this nucleon 
located at the intersection of the "external" Coulomb tra- 
jectory (corresponding asymptotically to the specific im- 
pact parameter) with the "working sphere" (see above). 
The same procedure is used to connect the direction of an 
outgoing particle at the nuclear periphery and its asymp- 
totic direction. 

These three modifications can be considered as manda- 
tory. They do not introduce any fitting parameter. Val- 
ues of the parameters have been fixed once for all, largely 
inspired from known phenomenology. 



2. Emission of clusters 

An improvement of the INCL4.2 model, concerning 
this feature, had been already proposed in Ref. The 
implementation of this feature in INCL4.5 is somehow 
different, although the basic idea is the same: an outgo- 
ing nucleon crossing the nuclear periphery is supposed to 
be able to carry along other nucleons to form a cluster, 
provided the involved nucleons are lying sufficiently near 
each other in phase space. We first describe the present 
implementation in INCL4.5 and then comment upon the 
difference with the work of Ref. [l5[ . 

The features of the model for cluster production can 
be described as follows: 

1 . An outgoing nucleon arriving at the surface of the 
"working sphere" , whether or not it has made col- 
lisions earlier, is selected as a possible leading nu- 
cleon for cluster emission, provided its energy is 
larger than the threshold energy, otherwise it is re- 
flected. 

2. Potential clusters are then constructed. The lead- 
ing nucleon is drawn on its (straight) line of motion 
back to a radial distance 

D = Ro + h, (3) 

Rq being the half density radius, and clusters are 
built by searching nucleons which are sufficiently 
close in phase space (A's are excluded)^. Clusters 
of increasing sizes are built successively. All po- 
tential clusters up to a maximum size A™"^ are 



^ If the line of motion of the nucleon does not cross the sphere 
of radius i?o + the nucleon is moved back to the minimum 
distance of approach of the center of the nucleus. 



considered. The criterion of sufficient proximity is 
expressed with the help of Jacobian coordinates: 

ri^[i^i]Pi^[i-i] < ha{Aci), for i = 2,3, A^i (4) 

where and are the relative coordi- 

nates of i—th nucleon with respect to the subgroup 
constituted of the first [i — 1] nucleons (i = 1 cor- 
responding to the leading nucleon) and where Ad 
is the mass number of the cluster. The value of 
hQ{Aci) is discussed below. The test on Jacobian 
coordinates is preferred to the usual test on the 
relative coordinates r^, pij for any pair of 
particles, because it disfavours exotic shapes (such 
as spaghetti) of the clusters. Considered clusters 
up to Ad""" =12 are listed in Fig. [H For the mo- 
ment, due to extremely fast increase of the com- 
binatorics with the mass of the cluster and due 
to limitations in computing time, clusters up to 
^maa: considered. Extending A^"-"" beyond 

this value may change the yields. We have checked 
that problems of convergence of the procedure man- 
ifest themselves by a slight excess of the yield for 
clusters with a mass number equal to as in- 

dicated in Ref. [13] , so that in practice numerically 
stable results are obtained for clusters of mass up 
to 7, and that the predicted yield for production of 
mass 8 can be considered as an upper bound. 

3. The less "virtual" cluster is selected. Let ^/s be the 
cm. energy of the composing nucleons, built on the 
4-momentum of a cluster, defined as the sum of the 
4-momenta of the nucleons inside the cluster. Let 
us consider the quantity 

= Xl"^')/^^' (5) 

where Bd is the (nominal) binding energy of the 
cluster. The cluster with the minimum value of 
V is selected. This quantity can be viewed as the 
excitation energy per nucleon of the cluster dimin- 
ished by twice the binding energy per nucleon. The 
introduction of this quantity is largely phenomeno- 
logical and is solely justified by the relative success 
of the model. 

4. The selected cluster is emitted provided three con- 
ditions are satisfied. First, it should have sufficient 
energy to escape, i.e. = ^(7i - Vi) - Bd >0 
, where the Ti 's are the kinetic energies of the nu- 
cleons and where the Vi 's are the depths of their 
potential wells. Second, the cluster has also to suc- 
ceed the test for penetration through the Coulomb 
barrier. Third, the cluster cannot be emitted too 
tangentially. If 9 is the angle between the direc- 
tion of the cluster (defined as the direction of its 
total 3-momentum) and the radial outward direc- 
tion passing by the center of mass of the potential 
cluster, it is required that 

cos e > 0.7 (6) 
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The idea behind this condition is that when a 
nascent cluster spends too much time in the nu- 
clear surface, it likely gets dissolved. These choices 
are admittedly made to improve the results at low 
energy, though some supporting arguments can be 
produced. 

5. If these tests are successful, the cluster is emitted 
with the kinetic energy Td in the direction of the 
total momentum of its components. If they are 
not, the leading nucleon is emitted alone provided 
it succeeds the test for penetration of the Coulomb 
barrier. If not, the leading nucleon is simply re- 
flected. 

6. At the end of the cascade process, short-lived clus- 
ters with a lifetime less than 1 ms (e.g. '^Li) are 
forced to decay, isotropically in their cm. frame. 
Clusters with a lifetime larger than 1 ms are consid- 
ered as detectable as such, prior to decay. Details 
are summarized in Fig. [TJ 
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FIG. 1: In this graph, the considered clusters for A^°'^ — 12 are 
identified by their charge Z (vertical ordering) and neutron 
A'' (horizontal ordering) numbersand are displayed inside the 
perimeter delineated by the heavy line. Absolutely stable 
clusters are denoted by the ordinary symbols inside yeloow 
cells. Cells displaying a time correspond to clusters with a 
lifetime larger than 1 ms, that are considered as detectable 
clusters. Other cells correspond to clusters with a lifetime 
smaller than 1 ms. They are forced to decay by the indicated 
mode. 

In our previous work [l5| , which was primarily devoted 
to production of clusters (up to Acz=4) at high energy 
(i.e. from 500 to 2500 MeV), two parameters, h (Eq. [3) 
and a single parameter Hq for all clusters in Eq. |4l with 
the selection of the heaviest possible cluster, were suffi- 
cient. In order to have a more or less satisfactory descrip- 
tion at low energy, we have been forced to consider differ- 
ent values of the parameter Hq for light clusters {Ad < 4) 
and a simple law {ha — hi{Aci/5y^'^) for heavier clusters. 
In INCL4.5, the parameters have been determined once 
for all, by fitting data on a few illustrative cases. The 
values of the parameters are given in Table HI 

3. Modifications concerning the Pauli blocking 

A strict Pauli blocking is applied to the first colli- 
sion: the nucleons should lie outside the Fermi sea af- 
ter the first collision. In INCL4.2, the Pauli blocking is 



applied stochastically, according to the product of the 
final state blocking factors. Conjugated with the fact 
that constructing the target with nucleons at random 
generates events with a nonuniform Fermi sea (even if 
it is uniform on the average over events), this procedure 
has the drawback of allowing sometimes collisions which 
otherwise (i.e. in a perfectly filled Fermi sea) would be 
strictly forbidden. On the other hand, it allows to ac- 
count for surface effects and for effects of the depletion of 
the Ferrni sea as the cascade process evolves. It is found 
in Refs. 26', 35] that a good compromise is achieved when 
a strict Pauli blocking is adopted for the first collision and 
when the usual procedure is kept for the subsequent ones. 

4- Modifications for soft collisions 

A. Soft collisions. In INCL4.2, soft collisions (with c. 
m. energy \/s < ^/sq = 1925 MeV) are neglected. His- 
torically, this choice was made to avoid inconveniences 
linked with the raising NN cross sections at low en- 
ergy. More profoundly, the underlying argument states 
that soft collisions (with low momentum transfer) do not 
change significantly the energy-momentum flow in the 
system and that their effect is likely to be more reason- 
ably accounted for by the nuclear mean field. Further- 
more, changing the boundary between soft and hard colli- 
sions has no sharp effect. The reason is that, in spallation 
reactions, soft collisions occur mainly when the colliding 
nucleons are lying close to the Fermi energy and thus that 
these collisions are largely Pauli-blocked, as explained in 
Ref. [35| . This argument breaks down when a low en- 
ergy incident nucleon makes a collision in the nuclear 
periphery, where Pauli blocking is not very efficient. We 
thus decided to lower ^/sq, still trying to keep the results 
roughly equivalent and to save computation time. The 
new value is now ^^==1910 MeV. 

B. Special treatment on the first collision. For the first 
collision, we even lowered the value of y/sQ to the min- 
imum, i.e. twice the nucleon mass, in apparent contra- 
diction with the above arguments above. However, at 
low energy, only a few (1-3, on the average) collisions 
occur. Neglecting a soft collision, allowed however by 
Pauli blocking, especially the first one, may amount to 
neglecting the event. This may have dramatic effects on 
the total reaction cross section, since the latter involves 
all kinds of events, be them hard or soft. It should be 
stressed that this procedure does not change the results 
at high energy (say, above ^200 MeV). Indeed, in tis 
case, the first collision is always a hard one. Neglecting 
a subsequent soft collision is most of the time harmless. 



TABLE I: Values of the parameters of the cluster formation 
model. The parameters are described m the text: h is given 
in fm, ho and hi are given in fm MeV/c. 
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See Ref. [36j for more details. In addition, a special pro- 
cedure, named 'local E', is applied to the first collision. 
In INCL4.2, the momenta of the target nucleons are too 
large in the nuclear surface, near the turning points, since 
they experience constant square well potentials, instead 
of smoothly varying potentials (this is "the price to pay" 
to keep the simplicity of straight line motion in our code, 
see Ref. [ll| for details). In the 'local E' procedure, when 
two nucleons are selected for the first collision, their mo- 
menta are "corrected" , by replacing their values by those 
assumed by the nucleons in a smoothly varying potential, 
at the same positions. When testing a pair of nucleons 
for collision, the NN cross section is calculated with the 
corrected momenta. After the collision, if it occurs, mo- 
menta are "corrected back" to the INCL4.2 prescription. 
Once again, this procedure has no effect on the first colli- 
sion at high incident energy, nor when this first collision 
occurs in the bulk of the target nucleus. These two mod- 
ifications are instrumental to give the predictions of the 
total reaction cross section in agreement with experimen- 
tal data at low energy. 

5. Modification of the status of the participants 

If after a binary collision or after a A-decay, a nu- 
cleon (obviously a participant) has an energy smaller 
than the Fermi energy plus a small quantity ^, put arbi- 
trarily equal to 18 MeV as a first attempt, it is considered 
thereafter as a spectator (it can be "re-promoted" as a 
participant if it gets a suflaciently large energy transfer 
in a subsequent collision with a participant). This pro- 
cedure, which is inspired from the Isabel code 37-39], 
may be motivated by various considerations: in a nu- 
cleus, the Fermi sea is not sharply defined, nucleons can- 
not be localized with precision when their energy is low 
and correlations may render the difference between a nu- 
cleon above the Fermi level and a spectator rather fuzzy. 
In fact, there is no compelling argument in favour of this 
procedure. It is included here because it gives slightly 
better results in some cases, in particular for production 
of clusters at low incident energy. Unfortunately, it leads, 
in our case, to unphysical results in neutron spectra at 
very low incident energy, as it is clearly illustrated in 
Fig. [5] below. That is why this feature has been changed 
in the INCL4.6 version, see below. For convenience, we 
will refer to this procedure as the "back to spectator" 
recipe. 

D. Additional features included in the INCL4.6 
model 

1. Further modification of the status of the participants 

As just stated above, the modification described in Sec- 
tion III C 51 produces unpleasant results in neutron and 
proton energy spectra, which will be illustrated when dis- 
cussing comparison with experimental data. In INCL4.6, 



participant neutrons are considered as spectators whence 
their energy gets below the neutron emission threshold, 
corresponding roughly to ^ defined above equal to 7 MeV. 
The same procedure is adopted for participant protons 
when their energy falls below their emission threshold 
plus two thirds of their Coulomb barrier. The whole pro- 
cedure is certainly more acceptable, although it is not 
really justified on consistency arguments, but it is val- 
idated a posteriori by the results, as shown below (see 
Figs. E] and [H). 



2. Use of experimental thresholds for nucleon emission 

In our standard model, a participant nucleon (of type 
i — n,p) can be emitted when it hits the surface of the 
potential with a kinetic energy E larger than the poten- 
tial depth, or equivalently when the difference of energy 
E with the Fermi energy, E — Ep, is larger than the 
model separation energy Si . If it is emitted, the nucleon 
acquires an asymptotic kinetic energy E^o given by 

Eoo = E — Ep — Si (7) 

In our standard model, the values of Si are fixed for 
a given target nucleus and equal to the differences be- 
tween the potential depths and the Fermi energies, re- 
spectively'^. In INCL4.6, the energy E — Ep is compared 
to the physical separation energy Sf^^^ , taken from mass 
tables, for the emission from the actual nucleus, i.e. the 
target nucleus left over when the candidate particle is 
hitting the surface. Eq. [7] is replaced by 

E^=E-E^- Sfy' (8) 

The effect is expectedly small, except at low incident 
energy, where the precise value of the threshold energy 
for a reaction does matter, or when the evolution of the 
target in the course of the cascade reaches the border of 
available phase space. 



3. Modified value for Rmax 

Using a value given by Eq. [2] for the radius for the 
"working sphere" of Rmax — i?o + 8a is safe at high 
energy, since this value allows a sampling of impact pa- 
rameters concerning all reaction events: it encompasses 
the outskirts of the nuclear density by a value which is 

greater than Vint = \J C/^at/tt, that we denote for sim- 
plicity as the "range of interaction" , cr^Jv being the NN 



^ The definition of Si is slightly more involved for energy- 
dependent potentials. We simplify the presentation on this minor 
point. 
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total cross section at the incident energy per nucleon^. 
This is no longer true at low energy (from a few MeV 
to a few tens of MeV), where the "range of interaction" 
may become very large, because of the raising NN cross 
sections. We thus decided to use 

Rmax = Ro + 8a + nnt = i?i + 8a. (9) 

Accordingly, we have increased the maximum time 
of the cascade, which now corresponds to the time of 
passage of the incident particle through the "working 
sphere" along a diameter, when this time exceeds the 
usual stopping time, given by Eq. [1] 

^. Treatment of cluster-induced reactions 

In INCL4.2, an incident cluster (up to an alpha parti- 
cle) is considered as a collection of independent nucleons 
with internal Fermi motion superimposed to the motion 
of the incident cluster as a whole (see Ref. [lil), ad- 
justed in such a way that the sum of the total energies 
of the constituting nucleons is equal to the nominal en- 
ergy of the physical cluster. In other words, the cluster is 
replaced by independent on-shell nucleons with the cor- 
rect nominal total energy, but with an incorrect (smaller 
than nominal) total momentum. This approximation is 
justified at high energy, where the missing momentum is 
relatively small, and where this model gives more or less 
satisfactory results [ll|, EO, IHl- However, it is not re- 
ally appropriate for reactions at low incident energy. In 
this domain, several features should be taken properly. 
First the total energy-momentum content of the projec- 
tile should be preserved. Second, the collective motion 
of the projectile has to be respected at low energy and 
Fermi motion has to be progressively restored when the 
available energy allows it. Liberating Fermi motion at 
a too early stage (like in INCL4.2) may be harmful at 
low energy. Indeed a nucleon with a backward-oriented 
Fermi motion velocity larger than the collective velocity 
the projectile may fly at once in the backward direction, 
which is indeed unphysical. Third, emission of projec- 
tile spectators with the influence of their Fermi motion 
should be possible and dominated by geometrical prop- 
erties. This property, which is well established at high 
energy (43-i44i] ■ is expected to survive, to some extent, 
at low energy, at least for large impact parameters. Fi- 
nally, compound nucleus formation should dominate at 
low incident energy. To fulfill these requirements, we 
have implemented an ad hoc model, made of the follow- 
ing ingredients. 

a. Initialisation of the incident cluster. Nucleon mo- 
menta p'j and positions r',; inside the cluster are gener- 



^ In INCL, collisions are decided on a minimum of relative distance 
basis, which is equivalent to introducing an effective interaction 
with a range equal to ri„t ■ 



ated as before [ll| (note, however, that a special method 
is applied to ensure = and = 0)- At 

the beginning of the event, the cluster center of mass is 
positioned on the classical Coulomb trajectory in such a 
way that one of the nucleons is touching a sphere of ra- 
dius Rcoui- The latter represents the Coulomb barrier. 
The value of Rcoui is taken from the phenomenology of 
the Coulomb barrier heights and has been tabulated as 
function of the target mass for p, d, t, ^He and '^He pro- 
jectiles. It is given in the Appendix. Of course, for very 
large impact parameters and/or low incident energy, the 
Coulomb trajectory may miss the sphere of radius Rcoui- 
The incident cluster is then positioned at the minimum 
distance of approach. 

For the energy- momentum content, the following pro- 
cedure is adopted. Let e[ and p'j, the energy and the mo- 
mentum of the nucleons generated (as in INCL4.2) in the 

frame of the incident cluster (one has e[ = \l p\ + m|, 
^^p'j = 0). The nucleons are put off-shell by changing 
their energy into Ci = e'^ — v, defined by: 

^e. = ^(e^-t;)=M,„e, (10) 

where Mine is the exact projectile mass. When Lorentz- 
boosted with the incident velocity /3i„c — Pine /Mine, Pine 
being the nominal momentum of the projectile, the nu- 
cleons acquire 4-momenta {Ei,pi) such that 

Y,E, = W^ne.^V^ - Pne, (11) 

where Wine is the exact total energy of the projec- 
tile. This rather crude procedure aims at preserving the 
energy-momentum content of the cluster when replacing 
it by nucleons. The latter are then off-shell, which nat- 
urally accounts for their binding. Naively, the quantity 
V can be viewed as the potential necessary to bind the 
nucleons to the right binding energy of the cluster. At 
large incident cluster energy, the 4-momenta of the nu- 
cleons inside the projectile, are basically the same in the 
two schemes. 

h. Generation of the geometrical spectators. When po- 
sitioned as indicated above, the nucleons of the projectile 
are separated into geometrical participants and geometri- 
cal spectators, which are those nucleons whose direction 
of motion intercepts or not the working sphere, respec- 
tively For performing this separation, the direction of 
motion of the nucleons is assumed to be parallel to the 
collective velocity vector of the cluster. Spectator nucle- 
ons are put on shell with momentum pi, defined above. 



* This "geometrical spectator" term is used to distinguish these 
spectators, which are defined here by a geometrical criterion, 
from the usual definition of the spectators that are those nucle- 
ons, which, in a true reaction event, are intercepting the target, 
but avoid, by chance, collisions. 
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and energy E^^'^'^ ^ \l'Pi ^ ' ^^'^ frozen further 
on. In order to preserve the correct energy balance, the 
energy of the participants are decreased correlatively and 
equaUy. 

c. Treatment of the geometrical participants. The fol- 
lowing critera are used. If a participant has an energy 
lower than the bottom of the target potential well, the 
event is discarded. This happens very rarely when, due to 
fluctuations, the generated Fermi motion energy is very 
small, but occurs close to thresholds for emissions of par- 
ticles. 

If all the geometrical participant nucleons have a total 
energy above the bottom of the target potential and if 
one of these nucleons has a trajectory cutting the sphere 
of radius defined in Eq. [S^nd if one of these nucleons 
at least has an energy below the Fermi energy, a com- 
pound nucleus is formed with all the geometrical par- 
ticipants and the target, with an energy which is equal 
to the available energy. The so constructed compound 
nucleus is ready for de-excitation (no further cascade is 
performed). Of course, its energy should be larger than 
the nominal ground state energy, otherwise the event is 
discarded. This choice is inspired by the fact that a com- 
pound nucleus is expectedly formed at low energy when 
the inhibiting effects of the Coulomb barrier are over- 
come. 

Finally, when all the geometrical participants have an 
energy above the Fermi energy, the usual cascade is ap- 
plied. However, before colliding for the first time, geo- 
metrical participants are propagated with the velocity of 
the incident cluster as a whole. Right before their first 
collision, they are given back their Fermi motion and are 
put on shell. We take Fermi motion into account for 
calculating the cross section and the kinematics of the 
collision. 

d. Coulomb polarisation of incident deuterons. When 
an incident deuteron is selected, it is positioned initially 
with one of its nucleons touching the sphere of radius 
Rcoui- For heavy targets, if this nucleon is a proton, it is 
interchanged with the neutron. This procedure mocks up 
the polarisation of deuteron by the Coulomb field of the 
target. It is often advocated that this effect is necessary 
to describe properly deuteron- induced reactions [45| . 

The model described above may appear as a rather ad 
hoc procedure. It has been inspired by our willingness 
of coping with several aspects of the low-energy dynam- 
ics. But it can compete with other models, like fusion 
models. It manages Coulomb effects and Fermi block- 
ing effects at low energy, but in addition, it leads auto- 
matically and dynamically to departures from compound 
nucleus formation, as the incident energy increases. Of 
course, we are well aware of the fact that a valid justifica- 
tion of our procedure can mainly come from a successful 
confrontation with experimental data. 

We have presented here the initializations of the 
nucleon-induced and cluster-induced events separately. 
We want to stress that these preparations are not disjoint 
however. Actually, the initialisation of a cluster-induced 



event reduces to the one of a nucleon-induced event when 
the cluster is reduced to a single nucleon without Fermi 
motion. 

III. COMPARISON WITH EXPERIMENTAL 
DATA FOR NUCLEON-INDUCED REACTIONS 

A. Introduction 

We will not make a thorough comparison with ex- 
perimental data, nor with the predictions of compet- 
ing models, such as Bertini [4^ l43 Isabel ISTI. [38| . 
CEM03 [is', '4^, BUU f505, IQMDfsT], JQMD [52] and 
the BIC cascade of Gcant4 [4] (see Ref . 53] , for a recent 
description of most of these models), since such a com- 
parison can be found in the IAEA Intercomparison [l^ 
for INCL4.5. We will restrict ourselves to a comparison 
with data in cases where the influence of the modifica- 
tions brought by INCL4.5 and INCL4.6 are important 
or when the physical significance of the modifications 
brought by INCL4.6 in comparison with INCL4.2 can 
be tested. For nucleon-induced reactions, the differences 
between INCL4.5 and INCL4.6 predictions are rather of 
minor importance, except for low-energy neutron spec- 
tra, which is illustrated in Section IIIIC2I AU INCL4.5 
and INCL4.6 predictions have been obtained with the 
adjunction of the ABLA07 de-excitation model (actually 
the ABLA07V5 version), described in Ref. ^54^. 

B. Total reaction cross sections 

This observable is entirely determined by the cascade 
stage and even by the first collision in this stage. Total 
reaction cross sections calculated with INCL4.2, INCL4.5 
and INCL4.6 are displayed in Fig. [21 Whereas the pre- 
dictions are roughly the same for the three models above 
200 MeV, there is a dramatic improvement brought by 
INCL4.5 and especially INCL4.6 in the 10 MeV-100 MeV 
range. We recall that the diff'erence between INCL4.5 and 
INCL4.6 predictions is mainly coming from a too small 
radius of the working sphere Rmax in INCL4.5, which, 
especially at low energy, underestimates the number of 
interacting events for large impact parameters. It is re- 
markable that INCL4.6 is able to account for the bump 
of the cross sections appearing at a few tens of MeV and 
for its main properties, except for the Be target. The 
systematic parametrization of the data given in Ref. [s^l 
seems to indicate that reaction cross sections are not sen- 
sitive to structure effects and are smoothly varying with 
incident energy and target mass number, except for light 
nuclei at low incident energy, as it is illustrated by Fig. [51 

Above 200 MeV, where there is no real difference be- 
tween the predictions of the three models, the variation of 
the total reaction cross section with the incident energy 
closely follows the variation of the NN total cross section, 
a feature that comes naturally in Glauber models [H^ for 
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FIG. 2: Comparison of INCL4.2 (black), INCL4.5 (blue) and 
INCL4.6 (red) predictions with experimental data for the to- 
tal reaction cross section as a function of the incident kinetic 



energy T, for neutron-induced reactions on ^Be, "^^Fe and 
Pb (first row) and proton-induced reactions on ^Be, ^^C , 
56 and Pb (2d and 3rd rows). Data are from 



Refs. 



Fe, 



total reaction cross section. It is interesting to look at 
the origin of the differences between the INCL4.2 and 
INCL4.5 (or 4.6) predictions below 200 MeV. This is in- 
dicated by Fig. [31 where we have compared INCL4.2 pre- 
dictions with INCL4.6 predictions when one of the addi- 
tional features of this latter model is removed: Coulomb 
deflection, removal of soft collisions and the so-called 'lo- 
cal E' correction (see above). This Figure clearly shows 
that, in proton-induced reactions, the rise of the reac- 
tion cross section when the incident energy goes down 
from 200 MeV to 30-40 MeV, is due to the removal of 
the cut on soft collisions for the first collision, which in 
fact takes full account of the strong increase of the rais- 
ing NN cross section with decreasing incident energy in 
this range, when the 'local E' correction, defined above, 
is applied. This correction alone is already giving an 
important increase of the cross section, as indicated by 
the dotted curves in Fig. [31 Of course, at very low en- 
ergy, the reaction cross section for incident protons is 
decreased by the Coulomb repulsion, which leads to a 
global sharp decrease below the Coulomb barrier. In 
fact, it appears from Fig. [31 that it is important to take 
account of all the features which determine the probabil- 
ity of the first interaction. These effects are rather well 



identified in macroscopic models for the reaction cross 
sections (stI Issj . Similar considerations can be done for 
neutron-induced reactions. However, in this case, the 
fall-off of the reaction cross section at low energy is basi- 
cally due, in our model, to the Pauli blocking of collisions 
at very low incident energy. 
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FIG. 3: Comparison of INCL4.2 (black dotted lines) and 
INCL4.5 (full red lines) predictions with experimental data 
for the total reaction cross section as a function of the inci- 
dent kinetic energy T, in various systems. Dashed red lines 
correspond to INCL4.6 predictions when Coulomb deflection 
is neglected. Dotted (dot-and-dashed) lines correspond to 
INCL4.6 predictions when so-called 'local E' correction (soft 
collision suppression) are neglected. See text for detail. Data 
are from Refs. [55hf5?i |. 



C. Emission of light particles 

1. Neutron multiplicity in the cascade stage 

Neutron multiplicity is of crucial importance for nu- 
clear energy applications. For our purpose here, it is 
convenient, following Ref. [s^, to consider the multi- 
plicity of neutrons with a kinetic energy larger than 20 
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TABLE II: Neutron multiplicities (second column), obtained 
by integration of the experimental double differential cross- 
sections of Ref. f5U l in proton-induced reactions on ^'^^Pb nu- 
clei at three values of the incident kinettc energy T , compared 
with the predictions of various INCL-ABLA couplings. 





Exp 


INCL4.2 


INCL4.5 


TNCT.4 6 






ABLA 


ABLA07 


ABLA07 




T 


= 800 MeV 




2-20 MeV 


6.5 ± 0.7 


6.8 


6.33 


6.73 


20 MeV-E^a, 


1.9 ± 0.2 


2.5 


1.84 


1.91 




r -- 


= 1200 MeV 




2-20 MeV 


8.3 ± 0.8 


8.1 


7.8 


8.27 


20 MeV-E™,, 


2.7 ± 0.3 


3.1 


2.39 


2.48 




T -- 


= 1600 MeV 




2-20 MeV 


10.1 ± 1.0 


8.8 


8.6 


9.16 


20 MeV-E„,, 


3.4 ± 0.5 


3.7 


2.79 


2.90 



MeV, since this quantity is sensitive to the cascade stage 
only. Results are given in Table [TTl One can see that our 
predictions with the latest version INCL4.6 are globally 
satisfactory, although the agreement with experimental 
data is somehow marginal at 1600 MeV. The decrease of 
the E >20 MeV neutron muhiplicity from INCL4.2 to 
INCL4.6 is due to the "back to spectator" procedure. It 
can also be seen from Table HIl that the predicted neutron 
multiplicities in the 2-20 MeV range are rather well de- 
scribed. In this case, of course, the merit is to be shared 
by ABLA07 and 1NCL4.6, since the latter determines the 
excitation energy of the remnant, prior to evaporation. 




1 
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FIG. 4: Comparison of INCL4.2 (blue) and INCL4.6 (red) 
predictions with experimental data for neutron double dif- 
ferential cross sections in p+^^^^Pb collisions at 1200 MeV, 
as functions of the neutron kinetic energy T and at differ- 
ent angles, indicated on the figure. For the sake of clarity, 
the spectra at the various angles have been multiplied by the 
indicated constants. Data are from Ref. [66j . 



2. Neutron energy spectra 

These quantities arc also of importance for applica- 
tions. We will not comment very much on this point 
since extensive results can be found in the IAEA Inter- 
comparison and since globally, all our models are giving 
rather satisfactory results for neutrons with a kinetic en- 
ergy larger than 20 MeV, i.e. for those neutrons which 
are produced in the cascade stage. As an illustrative ex- 
ample, we show in Fig. |4]the predictions of INCL4.2 and 
1NCL4.6 for the neutron spectra in p+^'^^Pb collisions at 
1200 MeV. 

The INCL4.5 (not shown) and INCL4.6 models give 
essentially the same results, for all targets and all inci- 
dent energies, except below ~100 MeV (see below). The 
predictions of these models are somehow less good than 
those of INCL4.2. The difference particularly concerns 
neutrons with a kinetic energy between 20 to 50 MeV 
at angles between 10 and 50 degrees, where INCL4.5 
and INCL4.6 underestimate the cross sections. This defi- 
ciency partly comes from the implementation of the pro- 
duction of composite particles by coalescence (inexistent 
in INCL4.2), which in a sense "eats up" neutrons and 
protons in this energy range, and partly from the use of 
energy-dependent nucleon average potentials. Note that 



this deficiency is reduced at 800 MeV and seems to be 
slightly increased for lighter targets (Fe, Zr). It can also 
be seen from Fig. |4] that INCL4.6 is slightly better than 
1NCL4.2 for the high energy part of the spectra. 

We show in Fig. [5] the kind of results that we ob- 
tain around 60 MeV incident kinetic energy. The strik- 
ing feature is the "dip" in the neutron spectra around 
10 MeV produced by the INCL4.5 model. This orig- 
inates from the "back to spectator" recipe, explained 
in Section III C 51 This unphysical feature disappears in 
1NCL4.6, since the quantity ^ is then reduced to 7 MeV 
for neutrons and the same quantity plus two thirds of 
the Coulomb barrier height for protons. Dips are so re- 
moved from the nucleon spectra. Dips of this kind are 
not visible at high energy (say above 100 MeV, see for in- 
stance Fig.|4]) because they are hidden by the overwhelm- 
ing evaporation contribution. The theoretical shapes of 
the spectra above the evaporation peaks in Fig. [5] are 
rather satisfactory, although the degree of agreement is 
substantially smaller than at high energy (see Fig. \^ . As 
a matter of fact, the shapes are almost perfectly repro- 
duced at 800 MeV and are progressively departing from 
experiment when the incident energy is decreasing. 
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FIG. 5: Comparison of INCL4.5 (blue) and INCL4.6 (red) 
predictions for the neutron double differential cross sections in 
n+^^Fe collisions at 65 MeV (left panel, data from Ref. 
and in p+ Pb collisions at 63 MeV (right panel, data from 
Ref. [6i]). 



3. Proton energy spectra 

We spot some of our results in Fig. |6l Clearly, the kind 
of discontinuity observed in low energy proton spectra at 
62 MeV incident kinetic energy and due to the use of the 
"back-to-spectators" trick in INCL4.5, has disappeared 
in INCL4.6, as shown in the upper left panel. 

The INCL4.6 results are globally somehow satisfactory. 
However, systematic deficiencies can be detected. Cross 
sections are overpredicted for proton energy in the 10 to 
~ 50 MeV range at lowest angles (upper left panel), al- 
though results for heavier targets (not shown) are much 
better. This seems to come from an overestimate of 
the quasi-elastic component. Spectra at small angles 
and in the 40-100 MeV energy range are underpredicted 
at 175 MeV proton incident energy (upper right panel) 
and to a lesser extent at incident energy larger than 800 
MeV (lower right panel) . This defect seems to disappear 
around 800 MeV (lower left panel), where the agreement 
with experiment is remarkable. The origin of this discrep- 
ancy is not clear since it occurs in a region of the spec- 
trum neighbouring the quasi-elastic and quasi-inelastic 
peaks. The latter are believed to come from with a sin- 
gle NN collision, and thus from events with a low number 
of collisions. It is thus surprising that our model, be it 
INCL4.5 or INCL4.6, is able to describe reasonably well 
the quasi-elastic and quasi-inelastic peaks (at least above 
250 MeV) and not this few collision component, except 
around 800 MeV. 

Let us also mention that our calculations overestimate 
the proton yield in the evaporation region at high in- 
cident energy (see lower right panel of Fig. This 
may point to a deficient competition between the various 
particle emissions in the different versions of the ABLA 
model (or to an incorrect excitation energy of the target 
in INCL4.6). Let us however remind that the spectra of 







FIG. 6: Upper row: comparison of INCL4.5 (blue) and 
INCL4.6 (red) predictions with the experimental data (sym- 
bols) for the proton double differential cross sections in 
p+ Fe collisions at 62 MeV (left panel, data from Refs. [66 . 
l64i]) and in p+'^^Ni collisions at 175 MeV (right panel, data 
from Refs. \6^. [65i] ) . Lower row: same for the proton double 
differential cross sections in p+^"^Pb collisions at 800 MeV 
(data of Refs. [66 . [6?i] . left panel) and in p+^^'' Au collisions 
at 1200 MeV (data of Ref. right panel). 



evaporated neutrons, which provides the dominant chan- 
nel, in this energy range, are rather described by the 
ABLA model, as indicated by Fig. HI 



We finally comment on the comparison with INCL4.2 
(not shown here, see Ref. lUi]). The predictions of the lat- 
ter concerning proton spectra are slightly better at high 
energy, with a less pronounced deficiency at small angles 
and intermediate proton energy. At low incident energy, 
say below 150 MeV, the shape of the spectra are slightly 
better (see Ref. |26j), but the normalizations are notice- 
ably lower than those of experimental data, since the 
total reaction cross sections are badly underestimated, 
as noticed in Section fill Bl 



12 



TABLE III: Average multiplicity of the clusters produced in 
the cascade and m the evaporation stages, for two systems, as 
given by the INCL4.6-ABLA07 model. 



cluster 


p (1.2 GeV) + ^^''Au 


p (63 MeV) + ^°*P6 




case 


evap 


case 


evap 


d 


0.505 


0.166 


0.0335 


0.00035 


t 


0.124 


0.081 


0.0102 


0.00011 




0.046 


0.0067 


0.00147 


0. 


'He 


0.167 


0.352 


0.0141 


0.0136 



4- Light cluster energy spectra 

This is an important feature of INCL4.6, since INCL4.2 
can not accommodate emission of clusters during the cas- 
cade stage. We have made extensive calculations with 
our new INCL4.6 model for double differential cross sec- 
tions of d, t, '^He and a production and compared with 
representative experimental data of Refs [H, [6^ [6^ l68l - 
|73| covering a large incident energy range (62 MeV-2.5 
GeV) and a target mass range spanning from Al to Bi. 
Extensive results can be found in IAEA Intercomparison 
report [l^ . See also Ref. (T^I which contains results con- 
cerning excitation functions for production of He and H 
isotopes, using basically INCL4.5-I-ABLA. We just dis- 
play here some typical cases. 

We first present results for p-t-^^^Aw collisions at 1200 
MeV in Fig. [71 Globally, the agreement is fairly good, ex- 
cept for high energy deuterons at large angles. We remind 
that at this incident energy [l^, alpha particles are pre- 
dominantly produced in the evaporation stage, ^iJe's are 
dominantly produced in the cascade stage, deuterons and 
tritons are produced in both stages, but definitely more 
in the cascade, as indicated in Table Hill (see also Ref. 74 1 
for a discussion of these dominances and of their energy 
dependences). The impressive agreement between calcu- 
lations and data for a production is due to both INCL4.6 
and ABLA07 models (evaporation yield is determined 
by the excitation energy left after the cascade). On the 
other hand the good description of the high energy tails 
of all particle spectra is solely due to our cascade model. 
There is no practical difference between the predictions 
of INCL4.6 and those of INCL4.5. A comparison with 
INCL4.2 itself is impossible as cluster production is not 
included in this model. However, a similar cluster pro- 
duction model (up to alpha particles) was implemented 
in the intermediate version |15|, refered later as INCL4.2 
with clusters. The predictions of the latter are really sat- 
isfactory in the energy range (600-2500 MeV), in which 



it has been tested |15|, except for high energy a produc- 
tion, whose yield is underestimated by this intermediate 
version. Compared to INCL4.2 with clusters, INCL4.6 is 
significantly better on this latter point but the most im- 
portant improvement regards production at low incident 
energy and production of heavier clusters (see Section 
IIIIC 51 below). 







FIG. 7: Comparison of INCL4.5 (blue) and INCL4.6 (red) 
predictions for the double differential cross sections for 
deuteron (upper left), triton (upper right), '^He (lower left) 
and "^He (lower right) production in p^^^"^ An collisions at 
1200 MeV. Data from Ref. [661. 



Results for a slightly different system (p+^°^Pfe) at 
low energy are displayed in Fig. [51 At such an incident 
energy, the average excitation energy of the remnant is 
very small, evaporation of charged particles is strongly 
hindered and Icp's, even a particles, are predominantly 
produced in the cascade stage, as illustrated in Table 
mil So the comparison with the data of Ref. [s^l offers 
a severe test of INCL4.6 and especially of the cluster 
emission module. It can be said that, on the average, 
the model captures the magnitude and the shape of the 
spectra fairly well, although the agreement is not as good 
as at high energy. Of course, partial deficiencies can be 
seen. For d, t and ^iJe emission, cross sections are overes- 
timated above the expected value of the Coulomb barrier 
(~13 MeV for unit charge particles), and underestimated 
at energies close to the incident energy and small angles. 
The importance of this last deficiency can hardly be as- 
sessed, since this part of the particle spectra is partly 
coming from coherent processes, which are outside the 
scope of the cascade simulations. Finally, for a emis- 
sion, the height of the Coulomb barrier used in the cas- 
cade (^-^26 MeV) is definitely too large. It is larger than 
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— INCL4.6 p(63MeV^+-""l'l)^d+X 

— INCL4.6 Cascade ^5 




INCL4.6 p(63MeV)+-°VbVHe+X 
INCL4.6 Cascade 





5. Heavier cluster energy spectra 



INCL4.6 p(63MeVl+-"*PbVHe+X 
INCL4.6 Cascade 




FIG. 8: Comparison of INCL4.6 predictions (red) to the ex- 
perimental data (symbols) for the double differential cross 
sections for deuteron (upper left), triton (upper right), He 
(lower left) and "^He (lower right) production in p+^''^Pb col- 
lisions at 63 MeV. The blue curves give the cascade contri- 
bution (only the blue curve is visible when the blue and red 
results are identical). Data from Ref. [63i] . 



the Coulomb barrier height used in the ABLA07 ver- 
sion, around 15 MeV. The effect of this barrier is clearly 
seen in the a energy spectra shown in Fig. [5] because the 
evaporation of a particles, though not important, is not 
negligible. These considerations explain the presence of 
a dip in the spectra around 23 MeV. Such a dip is not 
present or not visible in the spectra of the other light 
clusters. 

We want to mention that the kind of agreement 
achieved in Fig. [5] crucially relies on criterion (Eq. [S]) . 
If this condition is removed, the yield of low energy clus- 
ters is badly overestimated. 

We cannot multiply Figures here, due to lack of space, 
but we can say that more or less the same kind of agree- 
ment with experimental data is also reached for target 
nuclei as light as '^'^Al. 




7Be iNCL4.(i 

.1 16deg INCL4.6 Cascade 
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FIG. 9: Double differential cross section for He (left panel) 
and ^ Be (right panel) cluster production in p+^^'^Au colli- 
sions at 1200 MeV. The predictions of INCL4.6 -I- ABLA07 
(red histograms) are compared with the experimental data of 
Ref. [731 (black points) at different angles, indicated in the 
panels. The green histograms indicate the cascade compo- 
nent. 



This paragraph refers to emission of clusters heavier 
than alpha particles. The separation between light clus- 
ters {Aci < 4) and heavier clusters may appear arbitrary. 
However, if the production of light clusters by a sort of co- 
alescence mechanism is more and more accepted, the one 
of heavier clusters by the same mechanism may appear 
doubtful. Here we want to point out some illustrative 
resuhs obtained with the INCL4.6 + ABLA07 model, 
including the dynamical coalescence model described in 
Section III C 21 We show in Fig. [9] the double differen- 
tial production cross sections of ^He and ^i?e clusters 
in p+^^'^Au collisions at 1.2 GeV. One can see that our 
model reproduces rather well the magnitude of the cross 
sections and the shape of the spectra. Notice that this is 
not an obvious result since, for instance, the ^He cross 
section is almost an order of magnitude larger than the 
^Be one. Fig. |9] also shows that these ions are predomi- 
nantly produced in the cascade, at least for kinetic ener- 
gies larger than ~20 MeV for ^He and ~30 McV for "^Be. 
The situation is different for the production of ^Li and 
"^Li, shown on Fig. 1101 These clusters are predominantly 
produced in the evaporation, but the large energy part of 
the spectra, say above 50 MeV, is explained by the pro- 
duction in the cascade. Notice that these results are less 
satisfactory than those displayed in Fig. |9l since our pre- 
dictions seem systematically too low around 50-60 MeV. 
We have calculated production cross sections for other 
targets and other energies. Because of lack of space, re- 
sults will not be given here and will be described in a 
futurepublication. Some partial results are contained in 
Ref. [34| . Let us just mention here that our model yields 
similar results as those displayed in Fig. [S) However, 
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at low incident energy (~200 MeV 17511 or less), most 
of the energy spectra are too hard [Tg]. This matter is 
in progress. Nevertheless, our present results are rather 
unique in explaining the production of heavy clusters by 
a dynamic coalescence model. At least, they give credit 
to the plausibility of this mechanism, even if the details 
may not totally correspond to our coalescence model of 

Section urea ' 



1(7.10 MeV)+l'l) 





FIG. 10: Same as Fig. [9] for ^Li and ^Li production cross 
sections. 



6. Pion energy spectra 

We will not expand very much on this point, since it 
is of small direct importance for applications (note how- 
ever that pion production involves A excitation, which 
influences sizably the energy flow at high incident en- 
ergy). We present some results in Fig.[TT] It can be seen 
that INCL4.6 is describing the data quite well, especially 
for negative pion production. These results present a real 
improvement compared to INCL4.2. The respective mag- 
nitudes of the TT^ and tt^ cross sections are well repro- 
duced, indicating a correct isospin dependence of the ele- 
mentary cross sections. The differences at low energy are 
coming rather from the introduction of the pion potential 
well and the Coulomb barrier effects. Whereas the pre- 
dictions are correct for tt^, an overprediction of the tt"*" 
cross sections appears slightly above the Coulomb bar- 
rier. This shortcoming was already present in INCL4.2, 
but is now significantly reduced. 



D. Residue production 

1. Introduction 

Residue production properties are determined by both 
the cascade and the evaporation stages, whose respec- 
tive influences cannot often be easily disentangled. Fur- 
thermore, the ABLA model which have been used with 




p{7.10 McV)+l'li 




FIG. 11: Comparison of INCL4.2 (dashed blue lines) and 
INCL4.6 (red lines) predictions with experimental data for 
inclusive n'^ (left panel) and vr" (right panel) double differ- 



ential cross sections in p+ 
are from Ref. 



'Pb collisions at 730 MeV. Data 



TABLE IV: Average value of the mass number A, charge num- 
ber Z , excitation energy E* (in MeV) and angular momentum 
J (m h units) of the remnant in p+^^^ Ph collisions at 1 GeV, 
as calculated by different versions of INCL. 



INCL version 


<A> 


<Z> 


<E* > 


< J > 


INCL4.2 
without clusters 


202.73 


80.51 


135.9 


16.81 


INCL4.2 

with clusters 


201.92 


80.13 


137.5 




INCL4.5 


203.0 


80.4 


166.4 


16.44 


INCL4.6 


202.7 


80.4 


157.5 


16.04 



INCL4.2 is not the same as the (new) ABLA07 model 
which have been used with INCL4.6 (and INCL4.5). This 
complicates the analysis. Therefore we will limit our- 
selves to a few significant results which may enlighten 
the merits and deficiencies of INCL4.6. More informa- 
tion can be found in Refs. I18l-l20ll. 



2. Residue mass and charge spectra 

We present in Fig. [T2]the predictions of INCL4.2 and 
INCL4.6 for mass spectrum with the experimental data 
concerning the illustrative case of p-f ^"*P5 at 1 GeV. The 
theoretical results have been improved on the low mass 
side of the so-called evaporation residue peak and have 
been deteriorated in the A=180-200 region and to a lesser 
extent in the fission peak. 

We now try, for this typical case, to disentangle the 
changes brought by the modifications of the our cascade 
model from those introduced by the changes in the ABLA 
model. In the changes from INCL4.2 and INCL4.6, listed 
in Sections III CI and III Dl the ones that are the most rel- 
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FIG. 12: Comparison of INCL4.2 (green lines), INCL4.5 (blue 
lines) and INCL4.6 (red lines) predictions with experimental 
data for residue mass spectrum in p+^^^Pfe collisions at 1 
GeV. Data are from Ref. W\. 



evant for the residue mass spectrum are the introduction 
of cluster emission and the "back-to-spectators" trick. 
They, of course, influence the mass, charge, excitation 
energy and angular momentum of the remnants. To 
enlighten the discussion we give in Table IIVI the aver- 
age value of these quantities for the different versions of 
INCL. It can be seen that the introduction of cluster pro- 
duction slightly increases the average excitation energy 
and slightly decreases the average mass of the remnant. 
The "back-to-spectators" trick produces, as expected, a 
strong increase of < i?* > (even when this trick is "cor- 
rected" as in INCL4.6) and a slight increase of the aver- 
age mass of the remnant. The increase of the excitation 
energy is likely at the origin of the better description of 
the mass spectrum at the low mass end of the so-called 
evaporation peak. It is presumably responsible for the 
underevaluation of the cross section in the A=185-200 
region. However, the competition with fission may play 
a role as well in this region. The various changes intro- 
duced in the ABLA code [13, llll precludes a detailed 
determination of the respective influences of the model. 
Roughly speaking, one can say that the fission proba- 
bility, for excitation energy and fissility parameter typ- 
ical of the system above, has decreased from ABLA to 
ABLA07. Further small changes in ABLA07 when cou- 
pled to INCL4.5 and INCL4.6 have increased the fission 
probability. Fig. [T^] seems to indicate that the fission 



probability is probably rather satisfactory and that the 
remaining lack of cross section in the A=185-200 region 
is due to other features. 

For the sake of illustration, we give the charge and mass 
spectra for other systems in Fig. [TS] One can see that the 
predictions of the INCL4.6-ABLA07 model are in a good 
global agreement with the experimental data. We com- 
ment on the remaining discrepancies. The so-called deep 
spallation products are understimated in p-t-^^Fe at 300 
MeV. This is also the case for the A'-^IS residues in the 
p+^^Fe at 1 GeV and for the A--160-170 residues in the 
p^_238j^ system at 1 GeV. This might indicate an under- 
estimation of the frequency of high excitation events. But 
this may as well indicate an underevaluation of the emis- 
sion of the large clusters in the de-excitation. At least 
this seems to be corroborated, for the p-|-^'^*C/ system, by 
the underpredicted yield of the A^20-40 residues, since 
this kind of residues are expected to be produced by evap- 
oration (beyond multifragmentation, which very likely is 
not operative in this system [STj). 

We have to underline the almost perfect description 
of the mass spectrum in p-f ^°^P6 collisions at 500 MeV. 
Only the A=207 yield is overestimated. 

Let us comment on the fission fragment distribution 
in the •p^'^'^^U system. Although the magnitude and the 
width of the fission peak are well reproduced, some fea- 
tures are not satisfactory. On the one hand, the shape 
of the fission peak is not so well accounted for, probably 
signalling inappropriate shell effects. On the other hand, 
the plateau in the A=200-225 region is somehow overes- 
timated, pointing probably to an imperfect evaporation- 
fission competition at large excitation energy^. 

Finally, let us notice that odd-even effects are too large 
in p-|-^^Fe, presumably due to a too large pairing used 
in ABLA07. 

As for the isotopic distributions, we restrict ourselves 
to show our results for the illustrative case of 
at 1 GeV in Fig. [141 The shapes of these distributions 
for fission products are very well reproduced by both 
INCL4.2-ABLA and INCL4.6-ABLA07 models, with a 
slightly better performance for the first model. The iso- 
topic distributions for evaporation residues (right panel) 
are well reproduced by the INCL4.2-ABLA model, except 
for the overall yield for the lightest ones. On the contrary, 
the INCL4.6-ABLA07 model does not suffer from this 
disease, but some deficiencies appear in the shapes of the 
distributions for the residues close to the target: the Ph 
isotopic distribution is underestimated around A=200, 
the Os to Hg {Z = Zt - 6 to Z = Zt - 2) distributions 
are overestimated on the neutron-rich side and the ^°^T^ 
(Z = Zt — 1,A = At — 1) production cross section is siz- 
ably overestimated (this was also the case in INCL4.2). 
These deficiencies seem to come from a too low excita- 



The dip around A=214 is due to very short-lived a-emitters 
which escape detection. 
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FIG. 13: Comparison of INCL4.2-ABLA (green lines) and INCL4.6-ABLA07 (red lines) predictions with experimental data for 
residue charge (1st and 3rd rows) and mass (2d and 4th rows) spectra in p+^^Fe at 300 MeV (upper left quarter), p+^'^^Pb at 
500 MeV (upper right quarter), p+'^ Fe at 1 GeV (lower left quarter) and p+ ?7 at 1 GeV (lower right quarter), respectively. 
Data are from Refs. [66l - t66i] . 



tion energy for those remnants that are located on the 
neutron-rich border of the so-called evaporation corridor 
(i. e. the region of the {N, Z) plane populated by the 
residues), close to the target. We have not been able 
to convincingly relate these deficiencies to the change of 
a specific parameter or to the removal of a specific hy- 
pothesis. Let us notice that the INCL4.6-ABLA07 model 
reproduces rather nicely the distribution for the Bi iso- 
topes {Z = Zt + 1) produced through {p,xn) reactions. 

Similar deficiencies occur in the other systems refer- 
enced in Fig. [T2] for isotopes close to the target, al- 



though these deficiencies have a smaller amplitudes in 
the p+'^^Fe system at 1 GeV and they are almost van- 
ishing in the same system at 300 MeV. 



3. Restdue recoil veloctty 

Longitudinal recoil velocity measurements have been 
performed in inverse kinematics experiments carried on 
at GSI. We present typical results in Fig. [TS] along with 
our predictions, for proton-induced reactions on ^^Fe and 
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FIG. 14: Comparison of INCL4.2-ABLA (green lines) and INCL4.6-ABLA07 (red lines) predictions with experimental data, 
for so-called isotropic distributions in p+^'^^Pb collisions at 1 GeV. The left panel shows typical isotopes produced as fission 
fragments and the right panel show typical evaporation residues. Data are from Ref. 



^°^Pb at 1 GeV. The predictions of INCL4.6-ABLA07 for 
the average values of the longitudinal recoil velocity of the 
residues agree quite well with the data, for mass losses 
AA = At - A up to ~50 for p+'^°^Pb and to ~15 for 
p+^^Fe (roughly one quarter of the target mass number 
for both systems). For the former system, the apparent 
lack of agreement for AA > 50 is presumably (or at least 
partly) due to large fluctuations in the calculation arising 
from very small cross sections (see Fig. fT^ . For the other 
system, the production cross sections remain sizable for 
very large mass losses (relative to the target mass), and 
Fig. llSl indicates that our predictions start to depart from 
experiment when mass loss gets larger than 20. 

The lower part of Fig. [15] indicate that the predictions 
of INCL4.6-ABLA07 concerning the fluctuations are in 
remarkable agreement with the experimental values, even 
for larger intervals of mass loss (predictions of the pre- 
vious INCL4.2-ABLA model are also remarkably close), 
for both systems. 

These results provide a severe test for the cascade 
model. The average velocities are largely determined 
by the cascade process and the subsequent evaporation, 
being basically isotropic (in the remnant frame), does 
not bring significant contribution. This is also true, to a 
lesser extent, for the r.m.s. fluctuations of the longitudi- 
nal velocity, except when evaporated particles are more 
numerous than the cascade ejectiles, i. e. when AA ^ 



40 for p+^'^'Pb and to -15 for p-f ^Ve. 

Our INCL4.6 cascade model is thus very successful in 
describing these observables and in picking up the basic 
features of the generation of the recoil. It is turning less 
good in describing the average velocities for large mass 
loss in p+^^Fe at 1 GeV presumably because this system 
is characterized by a strong depletion of the target, for 
which models like INC are becoming less suitable. 

For relatively small mass losses (^15 for Fe and <30 
for Pb), the contribution of evaporation to r.m.s. fluc- 
tuations remains small. In these conditions, which allow 
a test of the cascade alone, the average value of the lon- 
gitudinal recoil velocity and its variance are remarkably 
linear in the mass loss variable, both in INCL4.6 and in 
the experiment. As discussed in Ref. [ll|, this is an indi- 
cation of the recoil resulting from a diffusion process, ar- 
guably due to independent successive NN collisions. The 
quality of the predictions of our cascade model appears 
as a strong support of its basic premises. We will not 
elaborate on this point here. We will instead discuss a 
little bit the difference between the results obtained by 
INCL4.6 and INCL4.2. Definitely, the predictions of the 
former are improving on those of the latter for the av- 
erage velocity, whereas the differences are rather minute 
for the r.m.s. fluctuation. For the average velocity, we 
checked that more than half of the difference for the Pb 
case, and almost the whole difference for the Fe case. 
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FIG. 15: Recoil velocity of the residues. Comparison of 
INCL4.2-ABLA (green symbols) and INCL4.6-ABLA07 (red 
symbols) predictions with experimental data (blue circles), for 
the average velocity (upper row) and for the variance of the 
velocity distribution (lower row), as functions of the mass loss 
AA, equal to the target mass number minus the residue mass 
number. The left (right) panels refer to p+^^Fe (p+^^^Pb) 
collisions at 1 GeV. Data are from Refs. [?S.[6Qj. 



IV. RESULTS FOR CLUSTER-INDUCED 
REACTIONS 



A. Introduction 

As we said in the Introduction, INCL4.2 can accom- 
modate light clusters as projectiles, considering them as 
collections of on-shell nucleons with Fermi motion and a 
collective velocity adjusted in such a way that the total 
energy of the projectile nucleons (in the target frame) 
is equal to the nominal projectile incident energy. This 
method is justified at high incident energy (where the 
difference between the adjusted velocity and the nomi- 
nal velocity is small anyway). Indeed, in Refs. [Til [ssj. 
good results were obtained for deuteron-induced colli- 
sions around 1 GeV and departures from the simple addi- 
tivity of a proton and a neutron cascades were illustrated. 
Satisfactory results were also obtained for other clusters 
in Refs. [io, SII. At low energy, this simple method can- 
not be accurate. Even more, the adjusted velocity can- 
not be defined for incident kinetic energy lower than the 
binding energy of the incident cluster. To cure this situ- 
ation, we have introduced the method described in Sec- 
tion III D 41 Our interest is mainly motivated by appli- 
cations, in particular to accelerator-driven systems and 
spallation targets, as quoted in the Introduction. As an 
example, astatine isotopes can be produced in Pb — Bi 
thick targets through, for instance, (a, xn) reactions in- 
duced by secondary a-particles [13, ISa IS^l- Following 
this motivation, we will mainly devote our attention here 
to the residue production cross sections. 



seems to come from the introduction of the emission of 
clusters in the cascade. In some sense, this emission al- 
lows a larger mass loss (in the cascade) without changing 
very much the pattern of the collisions for energetic par- 
ticipants (which are supposed to contribute the most to 
the recoil), i.e. without changing the recoil velocity and 
the excitation energy. This leads for a given recoil ve- 
locity to a corresponding increased mass loss. In simple 
terms, this pushes the curve for the average velocity in 
Fig. [15] toward the right of the Figure. We could not 
isolate clearly the cause of the remaining difference. Pre- 
sumably it is due to the introduction of forbidding col- 
lisions below the Fermi level (item 3.B of Section III Cp 
and to the introduction of "the back-to-spectators" trick 
(item 1 of Section lll Pp . Both do not affect very much the 
collisions involving energetic participants, but slightly in- 
crease the excitation energy. So, they increase the mass 
loss for a given recoil velocity, pushing also the recoil 
velocity curve of INCL4.2 toward the right. 

For completeness, we mention that we made the same 
calculation for p+^^Fe collisions at 500 MeV (not shown 
here). The results are slightly better than at 1 GeV. 



B. Total reaction cross sections 



A set of total reaction cross sections is given in Fig.lTSl 
Since this quantity is slowly varying with target mass and 
even with projectile mass, for the same velocity, we mix 
data for neighbouring target nuclei. Although the ex- 
perimental data are rather scarce, one can see that our 
model roughly picks up the basic properties of the en- 
ergy dependence of the reaction cross sections, namely 
a rapid rise at low energy followed by a slowly varying 
plateau. The magnitude of the plateau value is par- 
ticularly well reproduced for incident a-particles. For 
the d-induced reactions, this plateau value is underesti- 
mated. This may be explained by the contribution of 
the Coulomb dissociation of the deuteron, which is not 
taken into account. Evaluations of this contribution per- 
formed in Refs. [H, l90l-[93| point toward a sizable cross 
section, of the order of 100-300 nib for the d — Pb sys- 
tem. This contribution is expectedly smaller for d — Fe 
and somehow negligible for a-induced reactions, because 
of the large binding energy and the compactness of this 
projectile. One may notice the spectacular improvement 
obtained by INCL4.6 especially for deuteron-nucleus re- 
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FIG. 16: Total reaction cross sections (in mb) for deuterons (upper row) and for a or He particles (lower row), on Ni and 
Fe targets (left panels) and on Pb and Bi targets (right panels), as functions of the incident kinetic energy (for deuterons) or 
of the incident kinetic energy per nucleon (for He and He). The predictions of the INCL4.2 (INCL4.6) model are given by 
the green (red) curves. In the last panel, the dotted curves refer to ^_ffe-induced reaction cross sections. Data are taken from 
Refs. jg&iga]. 



action cross sections. Of course, this is consistent with 
the better reproduction of the nucleon-nucleus cross sec- 
tion by INCL4.6 in comparison with INCL4.2, as disussed 
in Section IIII B[ due itself to a better treatment of the 
first collision at low incident energy. The predictions of 
INCL4.2 are closer to those of INCL4.6 for a-particles, 
compared to deuterons, because for peripheral collisions, 
the probability of interaction is increasing with the num- 
ber of nucleons of the projectile. Finally, the threshold 
behaviour is different in INCL4.2 and INCL4.6: in the 
former, it is mainly dictated by the value of binding en- 



ergy per nucleon of the projectile, wheras in the latter, it 
is mainly dominated by the introduction of the Coulomb 
deflection and of the fusion (compound nucleus) process. 

Let us add a comment on the plateau values. First 
of all, the maximum of the cross section is losely cor- 
related with the maximum which occurs in the proton- 
nucleus cross sections for the same target at low energy 
(see Fig. [2]). The value of the reaction cross section at the 
maximum for cluster-induced reactions is larger than the 
maximum value of the proton-nucleus cross section and 
reflects the size of the cluster, which allows interactions 
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for larger impact parameters. Let us finally indicate that 
the slow decrease of the cross sections in Fig. [1^] is pre- 
sumably related to the decrease of the NN cross section 
in this energy range. 



C. Residue production cross sections 
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FIG. 17: Bi{a,xn) cross sections for x—1, 2, 3 and 4, 
respectively, as functions of the a incident kinetic energy. 
The full black curves correspond to the predictions of the 
INCL4.6+ABLA07 model. The dashed red lines represent the 
results using the Bass fusion model, followed by the ABLA07 
model. Data are taken from Refs. [95i-f9?i]. 



As a first example of our results, we show in Fig. [TT] the 
cross sections for the production of At isotopes through 
^°^Bi{a,xn) reactions. Globally, the shapes and the 
magnitudes of the peaks of the cross sections are satisfac- 
torily reproduced. The correct position of the respective 
maxima indicates a good description of the competition 
between the successively opening channels in the evapo- 
ration. On the other hand, the good evaluation of their 
magnitudes indicates that the total reaction cross sec- 
tion is correctly given by the cascade model, as shown in 
Fig. [ini Fig. [T7] also shows the results of a calculation 
using the Bass fusion model ^H, coupled to the ABLA07 
de-excitation model. It can be seen that in these condi- 
tions, where fusion dominates, our model yields results 
which compare rather well with those of a very commonly 



used fusion model. It may be argued however that the 
Bass model is not really suitable for alpha particles. 

A little bit more detail concerning the ^°^i?i(a, n) re- 
action is give n in the first panel of Fig. [THl where data 
from Ref. |99| at higher energy have been added to those 
of Ref. [93. The high incident energy part of the ex- 
citation function is a good test of the cascade model. 
Indeed, we have divided the theoretical cross section in 
Fig. [18] into three components: the first one (in green) 
corresponding to events where the particles are ejected 
during the cascade stage^, the second one (in red) cor- 
responding to events with particles emitted during the 
de-excitation of the target remnant and the last one (in 
blue) for events in which the emitted particles are coming 
partly from the cascade and partly from the evaporation 
stages. The red curves can be viewed as giving the com- 
plete fusion (-l-evaporation) contribution whereas the two 
other curves can be interpreted as corresponding to in- 
complete fusion and/or pre-equilibrium processes. For 
the ^°^i3i(a,n) reaction, the mixed (blue) component is 
identically vanishing, since a single neutron is emitted 
either in the evaporation of the fused system or in the 
cascade stage. From the first panel of Fig. [THl it can be 
seen that the cross section above ~30 MeV incident en- 
ergy is well described by our model. The emission of a 
single neutron is then entirely due to the cascade mech- 
anism. The remnant, i.e. the compound nucleus, is then 
sufficiently excited to emit two neutrons, at least, with a 
large probability and thus does not significantly feed the 
In channel. The same is true for the Bi{a,p) reaction 
(shown in the central panel of Fig. [T5|l . In that case, the 
evaporation component is very much suppressed, even in 
the 20-30 MeV range, due to the Coulomb barrier. Fi- 
nally, the third panel of Fig. [TS| displays the situation for 
the ^°^i?i(a, np) reaction. Here, above ~70 MeV incident 
energy, the cross section is dominated by the emission of 
the neutron and the proton during the cascade stage. 
At lower incident kinetic energy, there is a very small 
component for the emission of the two nucleons in the 
evaporation stage, but at low incident kinetic energy the 
cross section is dominated by the mixed emission of the 
two particles. 

Another example is provided by Fig. 1191 which shows 
deuteron-induced reactions on ^°^Bi. The shapes of the 
cross sections are globally reproduced fairly well, though 
some details are missed. The decrease of the {d, 2n) cross 
section on the right of the peak is somehow too slow. The 
precise threshold behavior for {d,n) and {d,p) reactions 
is not exactly reproduced, in spite of the use of exact 
Q-values. Actually, the theoretical {d,n) and {d,p) cross 
sections are roughly the same. However, the experimen- 
tal (d, n) cross section is roughly 4 times smaller than the 
{d,p) cross section in the region of the maximum ('-^lO 
MeV). This discrepancy might be due to the Coulomb 
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Alpha tot kin energy (MeV) Alpha tot kin energy (MeV) Alpha tot kin energy (MeV) 

FIG. 18: ^°^Bi{a,n) (left), ^°^Bi{a,p) (center), ^°^Bi{a,np) (right) cross sections, as functions of the a incident kinetic energy. 
Each predicted cross section using the INCL4.6-ABLA07v5 model (black lines) is splitted into a cascade component (green 
lines), an evaporation component (red lines) and a mixed component (blue lines). See text for detail. The black lines give the 
sum of the three components and are visible when they are not identical to colored lines. Experimental data are taken from 
Refs. [96|, [9§|] for '^'^^Bi{a,n) and from Refs. i lOQ ] for the two other reactions. In the first panel, the data of Ref. [9§j are 
singled out by the large error bars. 



dissociation of the deuteron. We will discuss this point 
below, after the presentation of another case. 

Results concerning another heavy target are reported 
on in Figs. [201 and [H] The first one shows the various 
^^^Ta{a, xn) cross sections as functions of the a incident 
kinetic energy. One can see a remarkable agreement be- 
tween the predictions of our model and the experimental 
data. We also notice that the sum of the (a, xn) cross 
sections nicely exhausts the total reaction cross section. 

The splitting of each of the (a, xn) theoretical cross 
sections into three components as defined in Fig. [18] is 
displayed in Fig. [21] The bump of each cross section is 
mainly due to evaporation. For the (a, n) reaction, the 
decreasing part of the cross section above the bump is 
entirely due to cascade. There is, of course, no mixed 
component in this case. Starting from (a, 2n) and going 
to more and more emitted neutrons, the mixed compo- 
nent grows and become dominant above the bump region. 
This simply means that when the incident energy starts 
to correspond to this region, the excitation energy of the 
compound system becomes sufficient to emit a; -I- 1 neu- 
trons with a high probability, the (a, xn) cross section 
is sizably reduced and the underlying reaction mecha- 
nism resembles more and more to the standard spalla- 
tion reaction mechanism, involving a cascade stage, pos- 
sibly resulting in incomplete fusion, and followed by a 
de-excitation stage, both stages emitting particles. It is 
remarkable that both the shape and the magnitude of 
the various bumps and the trend of decreasing part of 
the cross sections are well described by our model. This 
means that our model is catching the main features of 
the compound nucleus formation at low energy and its 



progressive change into a mechanism involving so-called 
pre-equilibrium features. 

We show in Figs. [22] and [23] the results of our calcu- 
lation for "^i/e- induced reactions on ^^^Ta. If the trends 
are correctly accounted for, the agreement is much less 
satisfactory than for the ^^^Ta(a, xn) reactions. The 
cross sections for x=\ and x=2 are substantially over- 
estimated, and the one for x =3 is still too large by a 
factor 2. 

The ^^^Ta{^He,xn) cross sections are split in the var- 
ious contributions in Fig. [53] Several features have to 
be noticed: there is no evaporation component for x =1, 
the mixed component is increasing very sharply from the 
very opening of the cross section for x —2 and x =3. This 
component is dominating for x =2. 

Before interpreting our results, we want to comment 
on the shapes of the experimental cross sections. For 
{a, xn) reactions, both for ^^^Ta and ^^^Bi, the cross sec- 
tions show the typical pattern of raising and decreasing 
stages, replacing the {x — l)n cross section and then giv- 
ing place to the {x + l)n cross section. This is due to the 
positive threshold values (the opposite of the Q-value), 
increasing with x. For example, for the ^^^Bi{a,xn) re- 
actions, the threshold values are located at 15.56, 20.62, 
28.43 and 35.61 MeV, for x =1, 2, 3 and 4, respectively. 
On the other hand, for the ^^^Ta{^He, xn) reactions, the 
threshold values are -10.83, 3.88, 11.16 and 19.88 MeV. 
At the Coulomb barrier (roughly 15 MeV), the i^He.n), 
i^He, 2n) and (^He, 3n) channels are already open. That 
is why the cross sections, at least for a: =1 and 2, do not 
show a typical bump (as for '^'^^Bi{a,xn)), but rather 
a steady increase followed by a plateau. A compound 
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FIG. 19: Deuteron-induced reactions on Bi. The pre- 
dictions of the INCL4.6-ABLA07v5 model (fuU red hnes) 
are compared with experimental data (symbols) relative to 
^°^Bi(d,n) (squares), '^°^Bi{d,2n) (full dots) and ^°'-^Bi{d,p) 
(triangles) reactions. Cross sections are given as functions of 
the deuteron incident kinetic energy. Data from Refs. ^idh - 
il03| ]. 



nucleus formed at the Coulomb barrier has already an 
excitation energy of 25 MeV and has thus few chance 
to decay to the In exit channel and even to the 2n exit 
channel. 

We now come back to our numerical results for 
^^^Ta{^ He, xn). In accordance with the last remark, in 
our calculations (see Fig. [231 the i^He,xn) channel is 
fed by the cascade only. According to Section III D 4[ 
the three nucleons of ^He are lying above the Fermi en- 
ergy, at the beginning of the reaction. Very likely, the 
neutron is emitted freely or quasi freely (after a soft col- 
lision for instance), and the two protons are kept inside 
the target with a small excitation energy, otherwise a 
neutron would be evaporated. The (^He, 2n) channel is 
largely dominated by a mixed emission: presumably, one 
of the neutrons is emitted (almost) freely, whereas the 
two protons initiate a cascade process leading to a rem- 
nant sufficiently excited to emit a neutron. The {^He, 3n) 
channel presumably corresponds to a transition toward a 
sequence of emissions by evaporation. 

It seems to us that the overestimation in our model of 
the {^He,n) and (^ife, 2n), and perhaps (^ffe, 3n) cross 
sections indicates that the separation between compound 
nucleus and cascade regimes in our model, explained in 



FIG. 20: Ta{a, xn) cross sections for x=l (blue), 2 (green), 
3 (red) and 4 (black), respectively, as functions of the a inci- 
dent kinetic energy. The full lines correspond to the predic- 
tions of the INCL4.6-I-ABLA07 model. Data (symbols) are 
taken from Ref. fig^]. The dashed line gives the theoretical 
total reaction cross section. 



Section III D 41 is probaly too crude, somehow underesti- 
mating the fusion and overestimating the cascade cross 
section. 

Another possible interpretation calls for the Coulomb 
dissociation of the projectile. One can imagine that, at 
low energy, below or slightly above the Coulomb barrier, 
the '^He is dissociated, say in d + by the Coulomb 
field before it really hits the nucleus. If the deuteron es- 
capes, the accompanying proton "sees" a Coulomb bar- 
rier which is twice as small as the one seen by the ^He, 
but its kinetic energy is three times as small. The pro- 
ton is thus repelled. This phenomena would shift a 
part of the (^He,n) and {^He,2n) cross section to other 
channels. An a rgum ent in favour of this interpretation, 
quoted in Ref. [104| . comes from the fact that the sum 
of the experimental ^^^Tai^He,xn) cross sections does 
not exhaust the total reaction cross section. This effect, 
based on Coulomb dissociation, is expected to be roughly 
controlled by a parameter that can be losely defined as 
rj — {Zr)'^/B, where Z, r and B are the charge, the ra- 
dius and the binding energy of the cluster, i.e. increasing 
with the electric polarisability of the pro iecti l e an d in- 
versely proportional to its "stability" |35l Il05l4l07 |. On 
these grounds, one expects that the effects encountered 
for ^He are largely reduced for a-particles, as one may 
indeed observe above. However, the effect should be even 
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FIG. 21: Splitting of the calculated Ta{a,xn) cross sec- 
tions (black lines) for x—1, 2, 3 and 4, respectively, as func- 
tions of the a incident kinetic energy, into the various contri- 
butions: cascade (green), evaporation (red) and mixed (blue). 
See text for detail. Data (symbols) are taken from Ref. 



more pronounced for deuterons. But this does not seem 
to be the case for the example of Fig. [191 Of course, the 
effect is expectedly reduced for light targets. 

We now turn to illustrative results for light targets. 
Fig. [M] shows a comparison of our predictions with ex- 
periment concerning (a, n) reactions on ^^Fe and ^'^Ni 
targets. The agreement for the ^^Fe{a, n) cross section is 
quite good. There are two sets of data points for the reac- 



tion Nija , n). Our predictions are close to the data set 
of Ref. IllOj l . In addition, one can notice that our model 
is able to respect the difference between the cross sec- 
tions for the two targets, that are in fact rather similar to 
each other. The reason is that our model properly takes 
account of the different Coulomb penetrabilities and, es- 
pecially, of the real Q-values. The latter are -5.053 and 
-7.938 MeV, for ^^Fe{a,n) and ^°Ni{a,n) reactions, re- 
spectively. 

Another example is provided by Fig. [551 In contrast to 
the ^He +^^^ Ta case, iUustrated in Fig.lH the (^iJe, n) 
and {^He,2n) are now satisfactorily reproduced. Like 
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FIG. 22: Same as Fig.[20]for ^^'^Ta(^He,xn) reaction. Data 
(symbols) are taken from Ref. [i6^]. The dashed line gives 
the theoretical total reaction cross section. 



for the previous case, these reaction channels are open 
at the Coulomb barrier, though not "widely" open. The 
thresholds for the two reactions are 3.41 and 5.07 MeV, 
respectively, and the Coulomb barrier lies around 8 MeV. 
The fact that we reproduce (^He, n) and {^He, 2n) cross 
sections in this case and not in the ^^^Ta case is not in- 
consistent with the explanation in terms of Coulomb dis- 
sociation of the '^He in this latter case. Indeed, the prob- 
ability of this dissociation is expected to be very small for 
a ^^Co target. Fig. [55l also shows an i^He.p) cross sec- 
tion for ^^Cu (in view of lacking data for ^^Co a slightly 
different target is considered), for which our calculations 
reach a reasonable agreement with the data. 

The last case refers to the ^°Ni{d,n) reaction, illus- 
trated in Fig. [26l Although the opening of the channel is 
well described in our model, the cross section is notice- 
ably underestimated, in similarity with the ^'^^Bi case 
(see Fig. UM . This observation may be consistent with 
the explanation in terms of the Coulomb dissociation. 
However, this effect should in principle be much smaller 
^^Ni, as stated above. 

In summary, these examples and similar others that 
are not displayed here show that our model yields rea- 
sonable results for total reaction cross section and residue 
production cross sections at low energy. This observation 
seems to indicate that the model efficiently describes the 
probability of forming a compound nucleus and the pro- 
gressive appearence of pre-compound emission, despite 
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FIG. 23: Same asFig.EUfor ^**^Ta(^He,3;n) . Data (symbols) 
n <^ 41 



the crudeness of its ingredients. The model can even de- 
scribe threshold behaviour satisfactorily, owing to the use 
of experimental Q"Values. 



V. DISCUSSION. CONCLUSION 

We have presented here the updated and improved ver- 
sion (INCL4.6) of our standard INCL4.2 model. We re- 
call that the philosophy at the origin of this latter model 
leads to include as much known microscopic physics as 
possible, without relying on parameters. Of course, this 
model relies on assumptions, as reminded in the Intro- 
duction. 

The extension of INCL4.2 has been realized in two 
steps, giving birth to the INCL4.5 and INCL4.6 versions, 
described in Sections III CI and III Dl The results obtained 
by INCL4.5 being available owing to the Intercompari- 
son organized by the IAEA [l^, we have put here the 
emphasis on the results obtained by INCL4.6. We will 
not compare either with other INC models, since such 
a comparison has been done in Ref. [Tsj . concerning the 
INC (or QMD) models described in Refs. [13, [H EF 
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FIG. 24: ^^Fe{a,n) (upper panel) and ^'"Ni{a,n) (lower 
panel) cross sections. The full curves correspond to the pre- 
dictions of the INCL4.64-ABLA07 model and the symbols rep- 
resent the experimental data. Data are from Refs. [tOS - fH^] . 



l52l . IllTl . illSj l). It shows that our model, coupled to 
ABLA07, is one of the best combinations for the de- 
scription of spallation reactions, at least for the set of 
data proposed in the Intercomparison, basically nucleon- 
induced spallation reactions on key targets (maainly Fe 
and Pb) for incident energy spanning from ~60 MeV to 
3 GeV. See also Ref. 1201 fo r a short analysis of the Inter- 
comparison and Ref. |119l | for a similar intercomparison 
restricted to excitation functions of residue production 
cross sections. The attention has here been focused more 
on observables which are accessible by INCL4.6 only or 
which have been ignored by the Intercomparison, such as 
the recoil velocities and heavy cluster production. 

We first want to discuss the new features of our model. 
Part of the extension has been realized in the spirit 
of INCL4.2, namely by introducing known and well- 
established phenomenology. This bears on the introduc- 
tion of energy and isospin-dependent nucleon potentials, 
of pion potentials, of Coulomb deflection of charged par- 
ticles and of experimental Q-values, as discussed in Sec- 
tion III Dl The impact of these modifications on the re- 
sults is globally modest. However, they are important for 
some peculiar features. This holds for (p, tt) reactions [3l| 
and threshold behaviour of light cluster-induced reac- 
tions leading to the emission of a small number of nu- 
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FIG. 25: ^^Co(^He,n) (left), ^^CoCHe,2n) (center), ^^CoCHe,p) cross sections, as functions of the ^He incident kinetic 
energy. Comparison of the INCL4.6-ABLA07v5 predictions (red lines) with experimental data from Refs. ^i43i - fj-Mi] . 
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FIG. 26: ^'^Ni{d,n) cross section, as a function of the 
deuteron incident kinetic energy. Comparison of the predic- 
tions of the INCL4.6+ABLA07 model (red curve) with the 
experimental data of Ref. ^H-fii] . 



cleons (see Section HVCp . 

A second part of the extension, that we want to single 
out, consists in the introduction of the dynamical coales- 
cence model for production of clusters during the cascade. 
It is based on the simple idea that an outgoing nucleon 



can carry along other nucleons which are close enough in 
phase space. We insist on the fact that this idea departs 
from usual coalescence model by two features: it deals 
with phase space and not the momentum space, and it 
depends on the instantaneous distribution of the nucle- 
ons in the course of the cascade process and not on the 
final distribution. The importance of this difference is 
discussed in Ref. [l^. This model is performing rather 
well with a single proximity parameter (valid for all light 
clusters) at high energy, say above 500 MeV. The model 
had to be refined at lower energy, basically on two points: 
the proximity parameter had to vary, moderately, with 
the type of the cluster, and tangential emission of clusters 
had to suppressed. As far as the results are concerned, 
these modifications allow a real advance, since emission 
of light clusters (up to a-particles) are well reproduced 
by our model, down to 60 MeV incident energy. See 
Ref. [i3 and Section UrCl It is showed in Section Hinil 
that this model can be extended to emission of heavier 
clusters. We will comment on this point later on. 

Finally, the remaining part of the extension is more 
founded on recipes, admittedly to remove discrepancies 
with experimental data, than on solidly established phe- 
nomenology or pieces of models. It is accompanied by the 
introduction of ad hoc parameters or assumptions. The 
most intrincate aspect of this remaining part of the exten- 
sion is given by the new procedure for handling incident 
clusters. In INCL4.2, the incident cluster is considered 
as a collection of on-shell nucleons, gifted with a Fermi 
motion, and having a total energy equal to the incident 
total cluster energy. This approximation, which is pretty 
reasonable at high energy per nucleon of the incident 
cluster (compared to binding energy per nucleon and/or 
Fermi energy), is bound to fail at low energy. There- 
fore, in our opinion, it is necessary to account for the off- 
shellness of the nucleons while keeping total momentum 
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and energy to their nominal values. Another necessary 
feature appears to be the almost unavoidable compound 
nucleus formation at very low energy, at least for suffi- 
ciently small impact parameters. In addition, incomplete 
fusion and/or soft excitations of the projectile should be 
mad epossible at grazing impact parameters. The last 
ingredient of a satisfactory reaction model should allow 
the progressive growth of the pre-equilibrium processes. 

The surprising result of our study is that all these re- 
quirements can be met in a model using as basic quan- 
tities the 4-momenta of all nucleons and leading to re- 
sults in reasonable agreement with experiment. Even 
more, our model, whatever its ad hoc features, possesses 
the unique property of describing the gradual move from 
compound nucleus to pre-equilibrium processes with a 
unique set of assumptions. Usually, two distinct models 
are used, with sometimes an additional model to ponder 
the relative importance of compound nucleus and pre- 
equilibrium. Actually, we have achieved, in some sense, 
for cluster-induced reactions, the same objective as in 
our cascade model for nucleon-induced reactions, namely 
describing absorpion and so-called pre-equilibrium pro- 
cesses with the same cascade tools (propagation and 
binary collisions of individual nucleons). It should be 
stressed, however, that this objective has been realized 
on a rather empirical basis. 

Let us now turn to a discussion of the predictive power 
of our new model. We first first on nucleon, pion and 
residue production cross sections, either differential or 
global, for nucleon-induced reactions. These observables 
were already accessible by our previous model INCL4.2. 
We remind that some of the modifications included in 
INCL4.6 (even restricted to those which have an infiu- 
ence on the observables under consideration) may appear 
to have some ad hoc character, like the treatment of soft 
collisions or the "local E" trick, described in Section lll CI 
As it is often the case, when ad hoc or empirical modifi- 
cations are added to an already performing model, some 
of the results are improved and some others are slightly 
less good than before. For instance, the cascade part of 
the neutron spectra are slightly less good with INCL4.6 
than with INCL4.2, which gives an excellent reproduc- 
tion of the data (see Fig.lH). On the other side, INCL4.6 
generally gives more excitation energy than INCL4.2. As 
a consequence, the residue mass spectra have noticeably 
been improved (see Section IIII D 21) . A considerably im- 
portant improvement brought by INCL4.6 concerns the 
prediction of the total reaction cross sections, which has 
also allowed an improvement of the results at low incident 
energy, well below 200 MeV, the usually quoted limit of 
validity of the INC models. 

We now turn to the new potentialities of our improved 
model, which, in our opinion, outpace the improvements 
mentioned above and others discussed in the paper. We 
will shortly comment on each of them. 

1. Our improved model allows the production of clus- 
ters, owing to the implementation of a dynamical 
coalescence model. We have presented here (and 



in Ref. [T^) results for double differential produc- 
tion of light clusters, up to a particles, which are in 
rather good agreement with experimental data, for 
a wide range of incident energy and target mass. 
We have in addition showed that the production of 
heavier clusters with a kinetic energy larger than 
the typical evaporation energies can also be gener- 
ated by the same mechanism. In our simulations, 
which requires the production of all potential clus- 
ters, we are limited, by computational time, to A<8 
clusters. We have shown very promising results in 
Fig. ini but not all cases show a similar agreement. 
Nevertheless, we are convinced that our implemen- 
tation of the coalescence module is flexible and that 
a more involved search of the better parameters will 
give eventually an overall satisfying description. 

2. In Section IIII D 3[ devoted to the recoil of the 
residues, we have shown perhaps the nicest of our 
results. Recoil energies are small, but their knowl- 
edge is of uttermost importance for technological 
applications, since they determine the damage to 
materials under irradiation. We think it is an im- 
portant result, because the agreement with exper- 
iment is impressive, but also because the distribu- 
tion of the recoil velocity shows the fingerprint of a 
diffusion process, which is in concordance with the 
description of the reaction process as a succession 
of independent binary collisions. 

3. In Section lTVl we have set out our model for cluster- 
induced reactions. If at high incident energy the 
cluster can be viewed, in first approximation, as 
a collection of individual nucleons (with perhaps a 
procedure to handle spectators), it is crucial to in- 
troduce departures from this picture at low energy: 
off-shellness of the nucleons, geometrical spectators 
and participants, formation of compound nucleus 
and progressive appearance of pre-equilibrium pro- 
cesses. We have shown that all these features can 
be taken care of with individual nucleons and by us- 
ing criteria involving geometry and the 4-momenta 
of the nucleons. Even more, the results displayed 
in Section |IV] testify of the success of this proce- 
dure. Of course, rather global quantities have been 
considered and the model should be tested also on 
double-differential cross sections. But the results 
accumulated up to now are very encouraging. We 
want to draw the attention on the importance of 
this potentiality of the model for technical applica- 
tions. In a thick target bombarded by a high energy 
beam, many secondary reactions will be induced by 
clusters produced in primary interactions. 

We hope to have shown that the improvement of our 
INCL model, from INCL4.2 to INCL4.6, has generated 
a powerful tool for the description of nucleon and light- 
cluster induced reactions in a large domain of incident 
energy spanning from a few tens of MeV to 3 GeV. We 
stress that this description is based on a single micro- 
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scopic cascade model, even if some ad hoc prescriptions 
have been added, as we have discussed in detail. This 
addition has been made to improve the predictive power 
of the model in view of applications to thick spallation 
targets. Actually, the model is already included as an 
option in the Geant4 MCNPX and PHITS [mj 
transport codes. Nevertheless, we think that the model 
is still prone to further improvement. We mentioned in 
this paper several results which are not sufficiently satis- 
factory, like for instance the yield of isotopes close to the 
target, which have large cross sections and whose accu- 
rate predictions are important for applications. We are 
currently working on such topics. 



At < 10, Rcoui = Rmax (sBC Eq. [2]). The last equality 
also holds for incident protons and pions. 



TABLE V: Values of the parameters a, b, c. 



incident particle 


a 


b 


c 


d 


0.2565 


0.78 


2.5 


t 


0.2504 


0.58 


0.5 


'He 


0.5009 


1.16 


0.5 


^He 


0.5939 


1.64 


0.5 



VI. APPENDIX. PARAMETRIZATION OF Rcoui 



The quantity Rco 
is parametrized as 

Rco 



(in fm), defined in Section fll D 41 



I.UzZt 



aA. 



2/3 



(12) 



where z is the charge of the incident particle. The values 
of the quantities a, b and c are provided in Table fVl 
This is valid for target mass number At > 10. For 
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